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Abstract In order to tease apart proximate vs. ultimate sources of variation in reproductive strategy, studies have 
increasingly focused on populations rather than species as the unit of interest. The reproductive parameters of 
Phrynocephalus przewalskii (Agamidae) in different populations within the same phylogenetic clade were compared 
in this study. Female SVL, clutch size, egg volume and clutch volume varied significantly among populations. With 
increase in latitude, clutch size increased, while egg size decreased. Relatively fewer but larger eggs were produced 
with increasing of population density. Food availability had positive effects on clutch size, but no effect on egg size. 
Our result indicated that latitude, food availability and population density may be the proximate factors affecting the 


reproductive parameters of P. przewalskii. 
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1. Introduction 


A key objective in the research of life history evolution 
is to understand the ecological factors favoring different 
life history strategies (Partridge and Harvey, 1988). Life 
history strategies often exhibit considerable variation 
at different levels: among individuals within the same 
population (e. g., Ballinger, 1977; Schwarzkopf, 1992; 
Abell, 1999), among populations of a single species 
(e. g., Du et al., 2005; Jin and Liu, 2007; Zuffi et al., 
2009), and among different species (e. g., Dunham, 1988; 
Zamora-Abrego et al., 2007). The variation in life history 
characteristics may arise from several sources, such as 
genetic divergence caused ultimately by natural selection 
and phenotypic plasticity induced proximately by 
environment factors. Teasing apart proximate vs. ultimate 
sources of the variation in reproductive parameters 
remains a major challenge for life history studies. 

As ectotherms, lizards highly depend upon ambient 
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climatic conditions and display substantial 
environmentally induced variation in life history 
features and, as such, have become ‘model organisms’ 
in evolutionary ecology (Vitt and Pianka, 1994; Pinaka, 
1999). An approach that has become popular in analyses 
of reptilian life-history evolution is applied and compared 
among conspecific populations (Ferguson and Talent, 
1993; Qualls et al., 1995). Previous studies suggested 
that the variation in life history strategies of lizards were 
the adaptive responses to various selective environments, 
such as altitude (Hancock et al., 1998; Jin and Liu, 2007), 
latitude (Wilson and Cooke, 2004; Du et al., 2005), food 
availability (Ballinger, 1977; Abell, 1999; Bonnet et al., 
2001) and climate (e. g., precipitation and temperature) 
(Adolph and Porter, 1996; Zuffi et al., 2009). 
Phrynocephalus przewalskii is a small toad-headed 
oviparous sand lizard which has a large distribution range 
in northern China and adjacent Mongolia (Urquhart et 
al., 2009). It has been extensively studied in systematic 
and ecological aspects (e. g., Xu and Yang, 1995; Liu et 
al., 1996; Urquhart et al., 2009). However, the patterns 
of variation in clutch size, egg size and mating behavior 
in different geographic and climatic zones, and the 
sources of the variation have not been investigated. The 
phylogenetic study on P. przewalskii with mitochondrial 
DNA revealed two major lineages, with clear population 
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genetic structure within both clades (Urquhart et al., 
2009). The clutch size, egg size, and mating behavior, etc. 
of nine populations from one group in Clade I (Urquhart 
et al., 2009) were compared in this study. Our main 
objectives are to provide detailed information of the 
variation in reproductive traits among populations, and 
to understand the relationship between environmental 
factors and reproductive features in P. przewalskii. 


2. Materials and Methods 


Specimens of P. przewalskii from nine sites around the 
Tengger Desert (Figure 1, Table 1) were sampled during 
the breeding season between April and June in 2008. 
Specimens were deposited at the Herbarium of Lanzhou 
University, were euthanized immediately and then fixed 
in 8% formalin. Before gonadal analyses, body mass 
and snout-vent length (SVL) of each specimen were 
measured. The body mass of females was weighed by 
electronic scales. The SVL of 261 reproductive females 
from the nine populations were measured by calipers. 
Clutch size (CS) was estimated from counts of enlarged 
vitellogenic follicles or oviductal eggs after dissection. A 
preliminary analysis showed that the estimated CS from 
counts of vitellogenic follicles and oviductal eggs had 
no difference (ANCOVA, P > 0.05 in all populations). 


Therefore, the data were pooled. Length and width of 
oviductal eggs were measured to the nearest 0.1 mm, 
and the volume was calculated with the formula for the 
volume of an ellipsoid: V=*/,/7ab’, where a is half of the 
shortest diameter and b is half of the longest diameter. 

Insect abundance was surveyed using pitfall traps 
(20 traps; 20x10x10 cm) placed in the habitat where the 
specimens were collected. The traps were set for a 24-h 
period and the mass of insect was weighed to 0.01 g. The 
population density of P. przewalskii was estimated by the 
number of individuals encountered during a field search 
which was done between 10:00 and 15:00 hours. Food 
availability was estimated by dividing insect mass by 
population density. 

All data were In or arc-sine transformed to achieve 
the conditions for using parametric tests. The mean value 
for SVL of females differed among populations and the 
mean clutch size varied significantly among populations, 
and both were analyzed by one-way ANOVA. The 
relationship of clutch size and body length (SVL) 
was evaluated by linear regression analysis. If linear 
correlations were not statistically significant, comparisons 
of inter-population characteristics were tested by one- 
way ANOVA. Otherwise, the one-way ANCOVA with 
maternal size as the covariate was used. The distribution 
of clutch size was different among populations and 


} y 


Tengver Desert 


| — ) 
Figure 1 Locations of the populations of P. przewalskii examined in this study. Shaded area represents the distribution of this lizard (Urquhart 


et al., 2009). Alax Zouqi (AZQ); Gantang (GT); Mingin (MQ); Shandan (SD); Alax Youqi (AYQ); Hexipu (HXP); Jingtai (JT); Dajing (DJ); 
Luanjingtan (LJT). 
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Table 1 Female reproductive characteristics of P. przewalskii in 9 populations within the same phylogenetic clade. All values are mean + 


number of observed individuals). 
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SE (n: 


. Altitude Latitude Population Food SVL . SVL Egg size 
Populations . od. Clutch size n 3 n 
(m) (°N) density availability (mm) (mm) (mm) 
AZQ 1448 38.85 30 0.231 52.9+0.6 3.1401 35 52.5£1.3 433.3417.7 8 
GT 1614 37.47 40 0.185 51.7+0.6 2.4+0.2 39 50.0+0.8 408.6423.7 13 
MQ 1369 38.58 19.2 0.223 50.3 + 0.6 2.4+0.1 38 497+1.3 44064205 12 
SD 1789 38.78 12 0.339 52.6 + 0.9 3.2+0.2 30 55.2£1,.3 352033.6 11 
AYQ 1491 39.22 12 0.589 51.7+0.8 3.3+0.2 34 52.141.1 32924162 14 
HXP 1704 38.39 5:7 1.644 50.0 + 0.8 3.0+0.1 29 52.1+0.8 427.6+28.9 10 
JT 1736 37.13 24 0.129 52.74 1.6 3.1403 16 55.442.6 425.8 +21.3 7 
DJ 1886 37.49 45 0.888 52.3 + 0.9 2.40.2 21 51.4+0.9 427.8+29.3 10 
LJT 1336 37.97 14 0.273 55.2+0.9 3.440.2 19 51.2+2.8 40044449 4 
Total 51.9+0.3 2.9+0.1 261 52.0+0.5 400.9+9.3 89 


was analyzed by Chi-square test. The correlation of 
clutch size with latitude and food availability, and the 
negative correlation with population density of different 
populations were analyzed by artial correlation analysis. 
Clutch size variation among populations was checked by 
Chi-square test. Partial correlation analysis was applied 
to check the correlation of clutch sizes with latitude 
and food availability, with sample sizes of different 
populations being controlled. Both clutch size and female 
SVL were In-transformed, and then were analyzed by 
linear regression to determine the effect of female body 
condition on the relative fecundity. Same statistical 
method was then used to determine the tradeoff between 
egg size and relative fecundity. All data were given as 


mean + SE and significance was set as P < 0.5. P < 0.05 
was considered as statistically significant. 


3. Results 


3.1 Female body size We measured 261 reproductive 
females from the nine populations (Table 1), and obtained 
the mean value of SVL of different populations (ANOVA, 
F252 = 2.97, P = 0.003). Pearson correlation analysis 
indicated that there was no significant correlation between 
female SVL and altitude (r = -0.024, P = 0.696) or 
latitude (r = -0.043, P = 0.490). 

Linear regression analysis suggested that female 
SVL had significant effect on the clutch size of eight 
populations (P < 0.022 in all case) except for Mingin (r = 
0.305, Fi 37 = 3.70, P = 0.062), while no significant effect 
on egg volume was found in any populations (P > 0.059 
in all case). When the data were pooled, clutch size was 
positively associated with female body size (r = 0.493, 
F; 259 = 83.30, P < 0.001), while egg volume was not (r = 
0.060, Fi 559 = 0.31, P = 0.578). 


3.2 Clutch size Clutch size of P. przewalskii ranged from 


1 to 6, and mean clutch size varied significantly among 
populations (ANOVA, Fy »;. = 6.37, P < 0.001, Table 1). 
When the effect of maternal body size was controlled 
by ANCOVA, the difference in clutch size among 
populations was still significant (ANCOVA, F; 251 = 6.78, 
P < 0.001). The distribution of clutch size was different 
among populations (Chi-square test, 73, = 70.529, df= 
32, P < 0.001, Figure 2). The prevalent clutch sizes of 
the populations, Alax Zouqi (AZQ), Shandan (SD), Alax 
Yougqi (AYQ), Hexipu (HXP), Jingtai (JT), Dajing (DJ) 
and Luanjingtan (LJT) were 2 to 4, while Mingin (MQ) 
and Gantang (GT) were 2 and 3 (Figure 2). 
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Figure 2 The distribution of clutch size 


Partial correlation analysis indicated that clutch size 
was positively correlated with latitude (r = 0.265, P < 
0.000, Figure 3) and food availability (r = 0.157, P = 
0.011), respectively, and negatively correlated with 
population density (r = -0.355, P < 0.001), but it was not 
correlated with altitude (r = -0.023, P = 0.708). 

89 female specimens which contained oviductal eggs 
were used in egg and clutch volume analysis. Significant 
population variation in the mean egg volume (ANOVA, 
Fy so = 3.19, P = 0.003) and clutch volume (ANCOVA, 
F; 79 = 2.05, P = 0.051) showed similar trend. Pearson 
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correlation analysis indicated that egg volume was 
positively correlated with population density (r = 0.217, 
P = 0.041), and negatively correlated with latitude (r = 
-0.276, P = 0.009, Figure 3), but it was not correlated 
with altitude (r = -0.028, P = 0.792) and food availability 
(r = 0.020, P = 0.854). Clutch volume showed a trend of 
negative correlation with population density (r = -0.206, 
P = 0.054), but it was not correlated with other factors (P > 
0.160 in all case). 


3.3 Trade-off between egg size and relative fecundity 
No significant relationship between egg volume and 
relative fecundity was found in any population (P > 
0.065 in all case). When the data were pooled, a trend of 
negative correlation was found between these two traits 
(r = -0.204, Fi s7 = 3.78, P = 0.055). The coefficient 
of the variation in clutch size (ranging from 24.5% to 
40.5%) was higher than that of egg volume (ranging from 
11.5% to 22.4%) for eight populations, except for SD 
population (27.6% vs. 31.7%). The overall coefficient of 
variance of the mean of the population clutch size was 
also larger than that of egg volume (32.4% vs. 21.8%). 
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Figure 3 Effect of latitude gradients on clutch size and egg size 
4. Discussion 


A significant positive correlation was found between 


female SVL and clutch size in P. przewalskii. Such a 
relationship has been revealed common in other lizard 
species with different clutch sizes (e. g., Du et al., 2005; 
Jin and Liu, 2007). But the female SVL could not explain 
all the variations in clutch size, because the difference 
was still significant after the female SVL was statistically 
controlled. The usual clutch size was different among 
populations, and the distribution of the clutch size varied 
significantly among populations. 

We failed to find any significant correlation between 
female SVL and egg volume in almost all the populations. 
The lack of correlation between female SVL and egg 
volume suggested that egg volume might be optimized 
by local adaptation. Our conclusion was further proven 
by the fact that the clutch size of P. przewalskii was more 
variable than its egg size which was in accordance with 
the prediction from the optimality models (which is based 
on the assumption that at any point in an organism’s life 
history there is an optimum percentage of available energy 
that should be diverted to reproduction to maximize the 
parents’ total contribution to future generations; Smith 
and Fretwell, 1974). Significant geographic variation of 
mean egg volume was found, which might be restricted to 
a specific population. 

High altitude characterized by short warm period, 
low mean temperature and great fluctuation in food 
availability (e. g., Kovshar, 1981) has profound effects on 
the life history strategies of lizards (e. g., Badyaev, 1997). 
Indeed, the species inhabiting high altitudes increase 
investment per offspring through producing small 
clutches of large eggs (Hancock et al., 1998; Jin and Liu, 
2007) in order to increase the survival of their offspring. 
However, no correlation was found between altitude and 
clutch size or egg volume in P. przewalskii. This might be 
due to the altitude gradient being too small to detect any 
differences. 

Latitudinal effect upon reproductive traits was 
another focused issue (e. g., Wilson and Cooke, 2004; 
Du et al., 2005). Clutch size usually increased with 
increasing latitude, while egg size decreased (e. g., 
Dunham, 1988). P. przewalskii inhabiting higher latitudes 
also produced relatively more but smaller eggs. Most 
hypotheses formulated to explain this pattern fall into 
two broad categories: the reproduction-survival tradeoff 
and the “lack” clutch hypotheses (e. g., Cooper et al., 
2005). The first is the ‘reproduction-survival tradeoff 
hypothesis’, in which females that allocate more energy 
toward survival must allocate less toward reproduction 
(e. g., Moreau, 1944). The second category is the ‘lack 
clutch hypothesis’, which assumes that clutch sizes are 
constrained by the amount of food parents could bring to 
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the brood (Ashmole, 1963; Cody, 1966; Royama, 1969). 
It is known that with the increasing latitude, predation 
decreases, while the overwinter mortality increases. 
Lizards residing at high latitude undergo considerable 
mortality in winter, and breed during the spring when food 
is also plentiful, while those inhabiting at low latitudes 
undergo heavier predation. In order to achieve high 
reproductive success, lizards produce relatively fewer but 
larger offsprings at low latitude, but relatively more but 
smaller offsprings at high latitude. 

Population density can be used as an indictor of the 
intensity of intraspecific competition to some extent. 
Under crowded conditions, lizards exhibit unique life 
history characteristics, such as small clutch size, delayed 
maturity and low frequency of egg production (Case, 
1983). It was known that larger offspring had a survival 
advantage under the conditions of increased predation 
and competition. When the intensity of intraspecific 
competition was high, P. przewalskii produced relatively 
fewer but larger eggs and increased investment per 
offspring by decreasing the total investment. Through this 
transformation of the reproductive strategy, the female’s 
total fitness might be enhanced by increasing the survival 
of offsprings. 

Low food availability for reptiles could result in the 
reduction in clutch size (e. g., Ballinger, 1977; Agrawal 
et al., 2001) and even delay in reproduction. We found an 
increased clutch size and clutch mass with increasing food 
availability, while egg volume did not vary significantly. 
The result from this study suggested that the females 
of P. przewalskii under good conditions would increase 
reproductive investment through increasing the number 
instead of the quality of their offspring. 

In summary, P. przewalskii clutch size increased with 
increasing of latitude, while egg size decreased. Relatively 
fewer but larger eggs were produced with increasing of 
population density. Food availability had positive effect 
on clutch size, but no effect on egg size. This study was 
restricted within the same phylogenetic clade, even within 
the same population group. So, we can say that latitude, 
food availability and population density may be the 
proximate factors which affect the reproductive characters 
of P. przewalskii. 
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